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ABSTRACT 

Cheng and Ho theory for the analysis of laminated compo- 
site shells has been used to study the action of combined pressure, 
axial load and torsional loads. The theory is more involved and 
gives improved results over Donnell's shell theory. 

The solution technique adopted does not allow the 
satisfaction of boundary conditions; though for long shells the 
end conditions do not appear to affect the buckling load to a 
large extent. The study is made for determining the buckling 
pattern by. solution of an eigenvalue problem. Specific values 
of m-number of half waves along axis of shell and n-number of 
waves around circumference of shell are found depending on 
whether loading is single or combined. The effect of laminate 
configuration, number of laminae and ply-orientations on buckling 
loads are studied for a shell subj ected to all combinations of 
loading. 

Optimization analysis is carried out for six-layered 
shells with ply- angle and ply- thickness ^s the design variables. 

The sequential augmented Lagrangian multiplier method is used 
for the constrained maximization of buckling load. Study is done 
for all cases of combined loadings with varying length-to- radius 
ratios of the shell. The effect of fixed lamina thickness and 
variable lamina thickness on maximum buckling load is also 
studied. It is observed that at the optimal point, the laminate 
configuration turns out to be very close to antisymmetric. 



Chapter 1 


INTRODUCTION AND LITERATURE SURVEY 


1*1 Brief History and Introduction : 

Shells have been known to be in use for.mankind since 
ages. The abundance of naturally occurring shell structures 
like bamboo, natural tunnel formations, bird's eggs etc. have 
indicated the remarkable strength properties and load carrying 
abilities of such structural elements. As a result of technolo- 
gically-inclined burgeoning of human interests, a purpose orien- 
ted study and experimentation on shell structures has created 
an area of scientific research which has continuously updated 
itself with time. 

The utility of shell structures like grain storage 
silos, tunnels, conduits, boilers, architectural domes etc. has 
highlighted the necessity of developing theories for accurate 
predictions of behavior of these structures. The parallel 
improvements in materials for use— and lately the composite 
materials — have led to highly enhanced applications of shell 
structures in previously unexplored realms. 

The period of 1920s saw the emergence Of experimentation 
in non-conventional structural materials, starting with light- 
weight, high strength alloys like duralumin, mylar etc. The 
potentials of such materials, and their drawbacks were realized. 
This put the designers on a new track of developing materials 
which would combine better, properties of the constituents. For 
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example, duralumin was found .fit by all structural strength 
demands, but failed where corrosion came into picture. This 
led to use of wooden chine strip over duralumin shell. It was 
the first crude attempt using a composite material, which 
ofcourse failed due to the lack of proper understanding of 
composite materials at that period in time. 

Aerospace structural design requirements of high 
strength-to-weight ratio and high stiffness- to-weight ratio, 
made the need for composite materials stronger. This acceler- 
ated experimentation on materials in companies, research labo- 
ratories and universities all over. World War II brought about 
an unprecedented upsurge in practical testing of existing theories. 
The ramifications in research in the field of shells went on as 
the quantum of knowledge leaped. 

Research on composite materials and composite struc- 
tures revealed the fact that their design had advantages of 
high strength-to-weight ratio, durability, ease of repair, design 
flexibility and high energy absorptivity. Also seamless cons- 
truction, as is needed for shells, was easily possible. With 
imperative demands calling up in the development of military 
aircrafts, jet engines, hypersonic flight and large-distance 
transportation, strides were taken towards research into 
chemistry and mechanics of composites. 

Research on the problem of buckling of laminated compo- 
site shells has been extensively done. With varying degree of 
approximation, the general theory has been improved upon. Every 
succeeding improvement is an attempt to reduce simplifying 



assumptions and thereby narrow the discrepancies with practical 
observations. 'Vith the availability of high speed computers, 
it is possible to use a theory employing smaller number of 
simplifying assumptions. Cheng and Ho [id theory is one such 
improvement on Donnell’s theory which has been extensively 
worked upon by many investigators. 

'I'S Literature Survey t 

Budiansky [lj has mathematically formulated the exact 
tensor equations of equilibrium of isotropic shells. Small 
perturbation theory to solve the non-linear shell equations has 
been used. This limits the validity of the theory only for 
small strains and rotations. 

Mah, Almroth and Pittner [2] have done a classical anal- 
ysis on the stability of cylindrical shells using Donnell’s 
assumptions. Ortho tropic as well as ring-stiffened shells have 
been considered. The loading applied to the shell has been any 
combination of axial compression, external pressure and torsion. 
The effects of these combinations on buckling loads of the 
shells have been studied. 

-Tennyson [3] has done. a photoelastic analysis of iso- 
tropic shells subjected to axial compression. This is done 
with study of high speed photographs (2000 frames/sec) of buck- 
ling process taken through a plane reflection type polariscope. 
Changes in isoclinic pattern with change in buckle wave shape 
have been recorded as a function of time. Predictions of 
theory, giving a buckling pattern of wave ^Lengths in 
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circumferential and axial directions have been verified by this 
photoelastic analysis. 

A small-deflection theory for general instability of 
ortho tropic circular cylindrical shells has been adopted by 
Becker and Gerard [4]. A solution for the governing differen- 
tial equations for chosen buckling load conditions has been 
shown to agree considerably with established data. Ring-stiff- 
ened cylinders have also been considered for the purpose of 
analysis and experimentation. 

Dong, Pister and Taylor [sj have postulated a theory 
for anisotropic plates and shells. They have taken the plate 
or shell to be composed of an arbitrary number of bonded layers, 
each possessing different thickness, orientation and elastic 
properties. Donnell-type equations have been discussed for 
determining stresses in individual laminae due to bending alone. 

Serapico [7] has investigated general instability of 
reinforced conical and cylindrical shells. The method of solu- 
tion is based on the principle of minimum potential energy 
applied to an equivalent orthotropic shell. General' results 
have been reduced to relatively simple formulae and are compared 
with other investigations and experimental values. 

Lowe [8] has presented a method for analysis of stresses 
in axi symmetrically loaded orthotropic filament wound cylinders. 
For a case of a circular cylindrical shell \vith semi-spherical 
dome, the stress analysis at dome-cylinder junction has been 
done. Numerical results have been compared to values obtained 
from hydro- test of a filament wound cylindrical shell. 
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Dow and Rosen [9] have evaluated the structural effi- 
ciency of various materials used in construction of orthotropic 
shells. A loading-index approach is used in which a measure of 
structural v,/eight is plotted as a function of an appropriate 
measure of design load. This is done such that the structure 
having the least value of ordinate at any value of abscissa is 
the one giving minimum weight for that design load. A non- 
dimensional structural efficiency parameter is obtained and 
plotted for different loading conditions. 

Weingarten [1 O] has studied the effects of pressure on 
small-deflection buckling of thin-walled cylindrical shells 
under bending by using modified Donnell's equations. Results 
indicate the variation of bending stresses with variation of 
pressure loading. 

Direct application of Donnell' s equations, adapted by 
Tsai [6] has been done by Khot [-11]. Principle of stationary 
value of potential energy is used for fomulating the buckling 
equation. Taking restrictions on radius- to-thickness ratio, 
variation of critical buckling load with configuration of ply- 
angles Is studied for a 3-ply laminated shell. He concludes 
that configurations of type -0, 0, 0 and -0, 7i/2, 9 at 0 hi 30 ^ 
give maximum critical buckling load and configuration of type 
0, -0, 0 and 9, -0, %/2 give minimum critical buckling load at 
0 RJ 900. The value of half wave length in circum.f eiential 
direction is found to be about 12-16 in all the cases. 

Cheng and Ho [l2]-[l3] used Fliigge' s equations [20] to 
develop a more accurate theory for buckling of anisotropic shells. 



Assuming a simplified single term solution, they have obtained 
•the solution and tried to satisfy four boundary conditions at 
the ends of the shell. By performing tests on ply-wood shells, 
the effectiveness of this theory has been pointed out. 

Jones and Morgan [14], using Donnell's approximations 
have obtained numerical results for cross-ply laminated circular 
cylindrical shells of boron -epoxy material. A relation bet- 
ween Batdrof shell curvature parameter given as 

.2 2 1/2 

z = ^ [{1 n - ^ } j 

^ ■ 

2 2 

and a buckling load parameter Kx = - N^L /% ^ is obtained. An 

investigation of Kx with number of layers of shell is done to 
make the results more practically applicable. 

Ugural and Cheng [is] have investigated the buckling of 
composite cylindrical shells under bending moment using the 
modified Cheng-Ho [12] theory. Coupling between in-plane stret- 
ching and bending is also studied. The numerical results give 
the variation of buckling load with length-to-radius ratio for 
the shell, and also the variation of buckling load with ply- 
angle orientation. Comparison with available data for plywood 
shells is made. The conclusion that is drawn is that under pure 
bending effect shells with face layers circumferentially orien- 
ted show more resistance to buckling than those face layers 
which are axially oriented. 

Hirano [ 18 ] has given a closed-form solution for 
axi symmetric and axially non-symmetric buckling of angle-ply 
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laminated circular cylindrical shells with many layers. Loading 
conditions, of uniforai axial and compressive loads have been 
considered. 

Whitney and Sun C1«J have used Cheng-Ho theory for 
arriving at an exact solution of buckling equations for an 
orthotropic laminated composite shell. It is found that though 
the series solution does not satisfy all the theoretical boun- 
dary conditions^ the experimental results with simply supported 
condition are comparable. 

Nshanian and Pappas [1?] have obtained optimum buckling 
loads for laminated composite cylinders. The governing equa- 
tions used are ones with Donnell' s approximations. Optimum 
buckling loads and corresponding optimum ply-angle orientations 
are obtained for cases of shells vdth different number of layers 
and different values of Batdrof parameter z. 

Hirano [I9j, using Donnell-Tsai approach, has considered 
the problem of optimization of cylindrical shells under axial 
compression. The shells are considered to be laminated with N 
orthotropic layers, and each layer having the same thickness and 
equal number of fibres in direction +a with respect to shell 
axis. By defining buckling stress as a function of buck ling- 
mode parameter the optimum results are obtained for a six- 
layered shell. 

^ of existing li terature reveals the foll- 

owing salient features: 

Use of simplified Donnell' s theory in solution of the 
buckling equations for shells provides less involved 
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solution, but at the cost of neglecting the terms repre- 
senting non-linear interaction between transverse shearing 
forces and rotations or twisting moments. 

^ A considerable discrepancy exists in the experimental 

results and numerical results obtained by Donnell' s theory. 

* Optimization is based on the simplified theories, with 
ply-angle kept as the only design variable. 

♦ Using an involved buckling theory, and a detailed solution 
technique, the comparative study between combined loading 
effects and other geometric and material properties has 
not been made so exhaustively. 

The following section gives a brief account of the 
work undertaken and a view of the development of the present 
thesis. 

1.3 Scope of Present Investigation : 

With the view of approaching practical results for 
buckling behavior of laminated composite shells, a detailed set 
of governing equations is used. These equations are solved 
using a series solution technique. The series terns are assumed 
such that for all combined loading conditions, simply supported 
boundary conditions are satisfied as closely as possible. This 
technique yields a large set of simultaneous equations depen- 
ding on the .terms of the series used. 

First part of the present work is the eigenvalue solu- 
tion of the set of equations, giving the lowest critical buck- 
ling load. The values of m (number of half-waves in axial 
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direction) and n (number of waves in circumferential direction) 
at which the lowest critical buckling load is obtained, give 
the nature of buckling pattern at those loading conditions. 

A direct comparison is made with the existing results 
and validity of our solution is ensured. During the course of 
this study it has been possible to obtain new results on comb- 
ined loading cases. Variation of critical buckling load with 
respect to angle of orientation of the plies and with respect 
to different number of layers is also studied. 

The second part of work is a study of optimization of 
six-layered Boron-Epoxy shells under the action of combined 
loadings. The design variables considered are fibre orientations 
and individual laminae thickness. Results are obtained for all 
combinations of loads for the following starting points: 

* The angles of orientation of each lamina as 0°, 30®, 45®, 

60® and 90®. 

♦ Cross-ply laminae. 

Symmetric angle-ply laminates with the angle configurations 
of [ 90 / 0 / 45 ]^,- [0/0/45]^, [30/45/90]^ and [l 5 / 45 / 6 o]^. 

The critical buckling load is maximized for preassigned 
geometry of the shell in terms of its length and radius. 

In Chapter 2, the development of shell theory which has 
been followed here is described. Chapter 3 deals with the 
solution technique and optimization technique employed. In 
Chapter 4, results obtained in the present work are presented 
in graphical form and discussed. Chapter 5 concludes the 
present study giving the major conclusions drawn. 



Chapter 2 


DHVSLOPMEMT OF SHELL THEORY 

A laminated cylindrical shell consists of a number of 
layers of thin orthotropic laminae, with the axes of symmetry of 
each lamina oriented at an arbitrary angle to the shell axis. 
Consider a laminated cylindrical shell (Figure 2.1), with its 
geometric axes oriented such that x-axis is defined along the 
length of the shell, z-axis is defined along the radius and 
0-axis is defined along the direction perpendicular to x and z 
axes. Origin of the axes is taken at the centre. The shell 
with radius ’a' and middle surface as the reference surface, is 
subjected to combined loading of external pressure p, axial 
compression P per unit of circumference and torque T per unit 
of circumferential length. 

2.1 Assumptions of the Theory ; 

Assumptions involved in the development of the theory 
are summarized as follows: 

1) The ratio of the thickness of the shell to the radii of 
curvature of its middle surface is small compared to 
unity. 

2) Displacements are very small as compared with the thick- 
ness of the shell. 

3) The straight fibres of an element which are perpendicular 
to the middle surface before deformation remain so after 
deformation and do not change their length. 
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Using these assumptions, the stress-strain relation is 
obtained from generalized Hooke’s law as: 
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( 2.1 ) 


where is the reduced stiffness matrix and the elements are 
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(2.2a) 
(2.2b) 
(2.2c) 
(2. 2d) 


In general, fibres in a lamina lie at an angle to the 
laminate reference axes. For the fibres aligned at an angle a 
to x-axis of reference plane (Figure 2.4), the transformed 
stiffness matrix for the lamina is 


=ij 


where. 
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The applicable Hooke's law now is - 
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From the equilibrium of the forces and moments acting 
on a face of an element of the K*” lamina (Figure 2.2), and 
defining z as the radial co-ordinate of the miodle surface of 
a lamina, the stress resultants and stress couples are given 

by Fliigge [20] as: 
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2.3 Strain-Displacement Relations : 

The deformation of the cylinder is described by the 
three components of displacement of any point on the midsurface 
of the shell. Defining u, v and w as the displacement of the 
midsurface in x, 9 and z directions respectively the strain 
relations are given by Cheng and Ho [l 2] as: 
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On expanding 1/(a + z) into power series as 

9 1 , 

1/zz^z'^ N 33 

a ^ - — + — - — . .) and neglecting z /a and higher 

B. 3. 

order terms, equations (2.7) become: 
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where the strains at the middle surface of the shell are: 
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= (u,Q/a) + V,, 


(2.9a) 

(2.9b) 

(2.9c) 


and the parameters corresponding to the change of curvatures 
are: 


X. 


= - w. 


xx 


Xq = - (l/a^)(w,QQ + w) 


^x0 


(2/a)[w,^Q + (u,Q/2a) - (v,y2)] 


0' 


(2.1 Oa) 
(2.10b) 
(2.10c) 


,th 


For K lamina of a shell, the stress- resultants and stress- 
couples are obtained from equations (2.6), by substituting equa- 
tions (2.8), (2.9) and (2.10) in equations (2.5) to get stresses 
in terms of displacements. For a laminated shell with N-layers, 
stress resultants and stress couples are obtained by summing 
for N laminae. Defining the following stiffness matrices: 
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N 


A. . 

13 

= S 

i=1 

(c. . ), 

13 'k 

~ \-1 ^ 

(2.lla) 

B. . 

13 

1 

" 2 

N 

.5 ‘=ij 
1=1 

^k ^^k " \-1 ^ 

(2.11b) 

°i3 

1 

“ 3 

N 

ifi ‘=13 


(2.11c) 


These are the following expressions: 


N. 


+ -4^) Ju,^} + + -P) { 


‘12 " a 


V, 


+ w } 

a a 


B., D., u,„ 


°ii V '^»ee ■*■ 

+ (B^, + "'xx’ ®12‘' 72 ’ 


^ (B +^) t ^ 5e^!lx, 

+ '®16 + 2a ' a 2 + a ' 

a 


(2.12a) 


N, 


9 




®12‘-'"'xx> ■" ^®22 - 


22 


^®26 2a a ^2 + a ^ 

a 


(2.12b) 


B 


N 


X0 


= ^6 + + ^A 


®26- -''■ 


— \ r ' 0 , W , 

26 ■*■ a ^ a a ^ 


, D, . u,Q 

+ ^^66 + ^’xi 
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D 


+ ‘®16 + + ®26'- 


w 


* 00 . ^ } 


^ ^ a5 


+ (B.. + 75^^) i - — T 

oo 2a a 


66^ . 119 + _!Z] 


2 a 


(2.12c) 


N, 


9x 


A„.?4fi + ¥)+ (A,^ + -a§)c-^ + v. 1 


"26 ' a • a' ■ • oo 


2a 


+ ■ “ )^- 


. '''^’90 , w 


W % 

2 ^ 

a Q 


+ (B 


66 2a 


D.. 2w, 

)C - — 


x9 ^ + Ilx } 


2 " a 


(2.12d) 


D. 


M 




11 a ' 'X 


'®16 + ■" ''•x' •" °1l‘-"’xx’ 




J ”'00 Si-) + D i- + _li) 

+ D^j!-— 2 - + 2 ’ + *^ 16 ' a ^2 a 

a ^ 


(2.12e) 


u. 


M 


0 


B)2i-'x> 




w 


»00 _ ^ 


22- ^2 -^32 


^»0 . ^»x 


^26 a " gS ^ a 


(2.1 2f ) 


M 


xQ 




■*■ ^®66 '*' a^* a ‘'■''’x*'*' °16*”*xxx ' 
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w 


+ DoaC 


’99 , w 


2w. 


26 ^ 2 2^ ^66 ^ " a 

a a 


x9 ’9 

5 

a 


U,Q V, 




(2.12g) 


K 


9x 




w, „ 

+ D,.{-w, ) + { ©9 + w } 

,16 ’XXX 26 2 ^2 

a ■ a 


2w. Q u. V, 


+ °65'-- a 


xQ X * X I 


2 ■ a 


(2.12h) 


2.4 Cheng- Ho Eguations 
Laminated ^ 


quations Governing Buckling of N-layered 
Composite Cylindrical Shell : 


The differential equations for shell buckling are given 
by Cheng and Ho [l 2] as: 


aN^ +.Nqv - P(^»QQ - ^’xQ = ° (2.13a) 


. ^ . ‘''Qv 

‘»x 


9 


XX 


x9 


^^9,q ^9,q ^^x9, " P®^^»QQ ■*■ '^’Q^ " a'^Pv, 


'9 


X 


2aT(v,^9 + w,x) = 0 


XX 

(2.1 3b) 


% + a(M^Q + + a^M^ - aNg - pa(au,^ - v,q + W,gQ) 

'ee ’XX 


- a Pw,^^ + 2a T(v,^ " ^.xG^ = ° 


(2.13c) 


Obtaining the derivatives of the terms N : 


% •’ %V 5 % ’ ^9 5 Mg ; Aig ; 

’9 ®^»0 >^®»x XX ’0 ^ ’ X ’99 x9 

“’'S-Qx 'denoting 

— pa / = k^q^ ^2 ~ ^/'^q2 ~ ^2^’ Pg ” "^^^22 ”” ^3^* 


22 ’ 
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the governing equations in terms of the displacements u, v and 
w for shell buckling are given as: 

2 

® “’xx ^'^11 ■*' ® 11 ^ “’xe ''15 “’ 6 S ^'^66 " ®66 ° 66 ' 

®^''-xx '^16 ®16 ° 16 ^ ® '''xo (^12 hb + ®12 + ® 66 > 


+ v.gg - a " '"-xx 6'='B,5 + D^g) 


™’xee^® 12 ''' ^65 " ^b() ~ '*'■ 600^*26 " °26 ^ 

+ "’ 0(^25 - ®26 t Dje' - '’'®^“-xx ‘^2 ^au.^gKj + 


K,(u,gg - aw,^)l = 0 


(2.14a) 


^ “’xx<h5 •" + D,^) + au,^„(5,, + 5^^ + + B,,) 


16 ^ “16 


x0' 12 66 66 ^ 12- 


+ u,qq(A 2 ^) + a + 2 D^^) + 2 av,^Q(A 2 . + 282 ^ 

+ D 26 > + ''■e 9 'l + B 22 > - ^®'*'-xxx <®16 + ^D,^) 

- + 2B, , - D_ + 3D^^) - aw.^ggCSB^^ + 2D^^ ) 


* xx9' 12 


66 12 


26 


" ""’©GO ^22 + + w. 


x' '26 " '■ 26' 


'9 


- q(a v ,^/2 + S'2av,^e + 2aw,^) + K,(v,gg + w.g)! = 0 


(2.14b) 


- ^^'^’xxx'^ll ° 1 l) - ^^“-xxe^^Bis + ° 16 ^ 

®'^’x 0 e ^®12 ■•■ ^^^66 ” ^ 66 ^ " “’ 900^®^26 " ^^ 26 ^ * ®“’x^^ 12 ^ 

+ u.qCSjj - ffjg + + 2lJ,|5) - a^v,^^g(B,2 
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°12 ^^66^ " ®^’X00^^®26 ''' ^^26 ^ 

- V, 999(622) + ^ ® 26 ^ ^’0 ^"^'''^’xxxx^^ll ^ 

+ "^^^'^'xxxO^^IS^ 2a2w,^^QQ(D^2 + 

+ x090^^26 ^ ''^»qq 90^^22^ ” ^ '^*xx^^®12^ 

" ^®'^’x0^^®25 " ^26^ “ '^’00^^^22 " ^°22^ 

+ w(l - B 22 + ^22^ ” 9iK^(v,Q - au,^ - ^, qq ) 

+ K2^®^'^»xx^ ^ ° (2.14c) 

2.5 Boundary Conditions; 


The governing equations (2.14) are complicated partial 
differential equations that are to be solved to determine the 
critical buckling load. The boundary conditions for shell with 
its edges supported tangentially and radially are given as: 


At X = + 1/2, 
w = 0 


u. 


B, , u, + D 

4 i 7 V 


1 1 ’ X 11 a 


B ^ 

X + _L2 V, 


© 


^12 ^12 ^2 
+ — v ,9 + w 

3 3 

B D , 

Sl6^»x - 7I ’^^0 

3 


+ °16'''x - 


XX 


£12 ^12 

' 2 '^*99 ~ _2 

3 . . , , . 3, ., , , 


E);, D. , D_ 

O J- 9 . w — 1^11 J. -la V - 0 

“ 2 a w»jje - 2 ^’Q a ^ 

3 , 


(2.15a) 
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Another pair of boundary conditions for the foregoing case is 
as follows: 


If V, u ^ 0 at both ends, then 

= 0 + TUjg/a = 0 at x = +1/2 (2. lob) 

If V, u = 0 at both ends, then 

v=0 u=0 atx= +1/2 (2.15c) 


If V ^ 0 and u = 0 at both ends, then 

+ Pv,^ =0 u=0 at X = +1/2 (2.15d) 

If V = 0 and u ^ 0 at both ends, then 

V = 0 N + lUj^/a =0 at x = +1/2 (2.15e) 

X y 

where = N^g + (M^g/a). 

2*6 Comparison v/ith Donnell's Theory : 


A major approximation according to Donnell's theory is 
neglecting all the terms in equations (2.10) except those invol- 
ving the second derivatives of w. This leads to a shortcoming 
that changes in curvature and twist are not fully accounted 
for. Further assuming the shell to be thin, the terms z/a in 
equations (2.6) are neglected compared to unity. Thus the ^ ^ 
expressions for strains of midsurface are simplified to: 



- z w. 


XX 


( 2. 1 6a ) 




u, 




w. 


@0 


W, 


- z 


X0 


( 2 . 16 b) 

( 2 . 16 c) 
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The stress- resultants and stress-couples are given as: 


N 

X 

h/2 
= / 
-h/2 

a dz 

X 

(2.17a) 


h/2 
= / 
-h/2 

CJg dz 

(2.17b) 

^X0 

= ^9x 

h/2 

= ^ Q dz 

-h/2 

-(2.17c) 

^x 

h/2 
= / 
-h/2 

z dz 

(2.17d) 


h/2 
= / 
-h/2 

Oq z dz 

(2.17e) 

''x9 

= '^^x. 

h/2 

= ^'vQ ^ 

-h/2 

(2.17f ) 


IVith these simplifications, the resulting governing equations are 
obtained on substitution of equations (2.16) and equations (2.17) 
in equations .(2.5 ). These equations are identical to the ones 
obtained by Jones [22j , and are mentioned here for a direct 
comparison with the modified form. 


® ■*' '^’89*66 ''' ® ''>xx*15 

+ av,^g(A,2 + Agg) + v,g^26 ' ^^'*'>xxx®11 
- - a'»',xee^®12 "*■ " *'>099^26 


+ '“•0<^26^ - “’xxS 


+ K,(u,gQ - aw,^)! 


(2.18a) 
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^ ^’xx^16 ®^’x9^^12 ^ ^ 60 '^ ^’99^^26^ '‘' 

'+ ^av.^gA^g + v,gg(l + Bjj) - - 

- - "'-9se®22 + “’xSe + “-Q 

+ K2(2av,^g + 2 aw,^) + K^(v,qq + w,q)} = 


a V, A , 
’ XX bo 




(2.18b) 


and 


^ ^’xxx®11 ^^^^»xx9®l6 " ®^’x99^®12 ^ “ '^»99Q^®26^ 


+ a'J»x'^12 ^»9^'^26^ " ^ ^’xxx^®l6^ " ^^^’xx9^®12 ^^66 ^ 


^^'^’x99®26 ^’999^22 "^>0 ^ 


X 26 


'9 


D. , 

xxxx 1 1 


■'■ "’xxxe^is ■*■ '">xx9e^'^i2 ■'■ ■'■ ‘*®"'’xeee^°26^ 


+ "'.9e0e°22 - 2a »'.xx®12 " ■*“’ xe®26 ' ^’ 99^22 * "'' ' " ®22^ 


2aw,^g)l = 0 


(2.18c) 


2.7 Note on the Present Theory ; 


Equations (2.14)are very general form of equations for 
buckling of a laminated cylindrical shell. In these, all the 
resultant forces and moments are taken into account, and any 
major assumptions which are likely to lead to deviation of the 
solution from actual performance, have been eliminated. From 
Donnell's equations (2.17a,b,c) it is seen that the term A^ ^ , 
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and appear separately, hence the simplification mentioned 
below is possible. In the equation 




— 

A * 3 

1 

r eO| 

(2.19) 


L w J 


B j D 

1 X J 


defining 

a = 

[Ar’ 



(2.20a) 


b = 

[aB j 



(2.20b) 

and 

d = 

[b^B 

- D] 


(2.20c) 

it is possible to 

obtain two equations of the form 




= [a] 

{N} + [b] { X } 

(2.21a) 

and 

{M} 

1 1 

II 

CN) + [d] CX } 

(2.21b) 


Solution of the equilibrium equations (2.13) using equation 
( 2 . 21 ) gives the critical buckling load for the shell. 

The- use of involved Cheng-Ho theory prevents such a 
solution procedure. Hence it is necessary to assume the displa- 
cements of the middle surface as some functions, which is discu- 
ssed in Chapter 3. 





FIG. 2.2 FORCE AND MOMENT EaUILIBRIUM FOR A SHELL ELEMENT 
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Chapter 3 

ANALYSIS AND OPTIMIZATION OF LA^/INAT5D S;-i£LLS 

A theoretical analysis for solving the buckling problem 
of a laminated composite circular cylindrical shell subjected 
to combined axial, radial and torsional loads is presented in 
this chapter. The governing equations (2.14) are based on the 
theory proposed by Cheng and Ho [l 2] for anisotropic cylindrical 
shells. By assuming a suitable displacement function for the 
midsurface displacements, the governing equations reduce to an 
eigenvalue problem consis.ting of a set of linear homogeneous 
algebraic equations. An effective numerical scheme is employed 
to perform the numerical com.putations for obtaining the lowest 
eigenvalue. 

As one would like to increase the load carrying capa- 
city of the laminated composite shell for the prescribed 
geometry and weight, optimization study is undertaken for six- 
layered shells subjected to axial, radial and torsional loads. 
The ply-angle and ply- thickness are kept as design variables 
and buckling load is maximized for the shells of prescribed 
geometry and weight. 

3*1 Selection of Displacement Function : 

The assumed displacement function for midsurface displa- 
cements should satisfy the governing differential equations 
(2.14) and boundary conditions (2.15). Since equations (2.14) 
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are complicated partial differential equations, it is impossible 
to arrive at an exact solution, using any form of displacement 
function. Hence a displacement function initially proposed by 
Flugge [20] and later tested by Cheng and Ho [12J, for laminated 
shell buckling is used in the present work. This displacement 
function adopted here can be used, for long shells, for all 
combinations of axial loads, external pressure and torsional 
loads. 


The midsurface displacements in x, 0 and z directions 
are assumed as: 


00 


u = E U sin(— + n0) 
. n a 

n=1 


00 


V = S V sin(-^ + n0) 
. n a 

n=1 


00 


cos(^^ + n0) 
n=1 


(3.1a) 

(3.1b) 

(3.1 c ) 


where 


X = 


mita 


(3.2) 


This solution describes a buckling mode with m. number of half- 
waves along axial direction and n number of waves around circum- 
ference of the shell. 

3.2 Satisfaction of Boundary Conditions : 

With the chosen displacement functions for u, v and w, 
the boundary conditions to be satisfied for a simply supporteo 
shell with origin at centre are given by equations (2.15). 
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For the condition w = M = 0 at x = +1/2, it is not 
possible to make all the coefficients of stiffness matrices 
Ai j , and zero for general ply laminates. The displa- 

cement conditions given by equation (2.15b) cannot be satisfied 
for combined loading conditions. However, if the length of the 
cylinder is very large, then the boundary conditions at the 
ends do not greatly affect the magnitude of critical stresses. 

Thus for the length 1 assumed to be a part of a much larger 
cylinder, the solution is not affected by end conditions. 

Ugural and Cheng [15] have solved equations (2.14) 
using a single series solution for the case of pure bending 
loads for 4-layered ply-wood shells and compared the results 
with experimental data. Similar procedure is followed here also. 

3.3 Solution Technique for Long Cylindrical Composite Shells : 

The governing equations (2.14) are solved for combined 
cases of loading, by substituting u, v and w in the form of 
equations '(3.1). For this it is essential to obtain the deri- 
vatives of u, V, w upto fourth order. Determining them and 
using in equations (2.14), one gets: 

S [ X.^ sin(~ + nQ) { A + B ^ ^ } + 2n X sin(~ + nQ) 
n=1 • 

+ n^ sin(^ + n9) - Sod +^6^ 7? sin(^ + n0) 

II^Y+ 2ff^^ + ^66^ ^12+ ^66^ 

+ sin(-^ + nQ) {A^,) + sin(^ + n0) } 
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+ A^n sin(M + ne) jjg + 0,^ 1 + An^ w sin(^ 


16 ‘ "16 ' •‘‘n 


^®12 - D..} + n^ sin(^ + n9) {] 


26 


"^26 ^ 


+ '^"'n =in(^ ^ n9) {a, 2 ? . „ 3l„(^ ^ nS) _ 3^, + ) j 


= q S U sin(^ + n9)K^ 

n=1 n a ^ 

+ 2An =An(~ + ne)Kj + {n^ + "9) 

+ AWr^ sln(^ + n0))K,j (3.33) 

00 2 

[A sin(A| + ne) {S,g + 2 B,g + + An sin(^ + n 0 ; 

**12 ^66 ®66 ®12* ■*■ ''®^ *^26* 

+ a 2 v„ sln( 2 f A n 9 ) + 35 ^^ + 20^,5 

+ 2 n X V^ sin(^ + n 9 ) } 

+ n^ V^ + n9) {1 + 3^2) + + ^8} 

^^16 °16^ - \ sin(^ + n9} {B.^^ + ^S6 “ ^12 ^^66 ^ 

"*" ^^®26 '*’ ^^26^ '*' ^n + n9) 

^®22^ + X sin(^ + n9) ^ ■*■ ^26^ " 

, ,, ;00 ' ^ ■ 

= q 2 [?^ V sin(-^ + n9) K^ 

n=1 ^ “ 2 

+ 2 X {n V sin(^ + n 9 ) + W sin(^ + n 93 }K„ 

+ n{nV^ sin(^ + n0)iKiJ 


(3.3b) 
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and 

cos(^ + nQ) + '^ 11 ^ + U^ cos(^ + n0) 

+ An2 cos(?^ + n©) f B ^ ^ 1 

+ Uj^ cos(^^ + n9) { B^^ - + X U^ cos(^ + nQ)' 

+ n cosO— + n9) 4 X^ v^ cos(^ 4 n©) 

{B . 4 2D .} 4 X^n V„ cos(^ 4 nG) C B, ^ 4 2B . . + D, ^ 4 3D, .. 1 
Id id n a 12 oo 12 oo 

+ Xn^ cos(^ 4 n0) ^33^^ 4 SD^^} 4 n^ V^ cos(~ 4 n9) ^^ 22 "^ 
4 X V^ cos(^ 4 n0) + B^^} + ^ 

4 x"^ W^ cos(^ 4 nG) {D^^} 4 4x^n W^ cos(^ 4 nG) fD^^} 

4 2 X^n^ cos(^ 4 nG) {D ^2 + 2D^6 ^ 

4 4Xn^V/ cos(^ + nG), {D_ } 4 n"^ ’// cos(^ + nO) {D_- J 
n a 2o n a zz 

4 x^ cos(^ + n9) ^28,, 2 ^ + 2Xn VV^ cos(^ 4 nQ) { 262 ^ - 
4 n^ \ cos(^ + nG) (2822 - 2022^ + cos(^ 4 nG) 

{ 1-822 + ^ 22 ? ] = q {X cos(-| 4 nG) 

+ n cos(^ 4 nG) 4 n^ cos(^ 4 n0)}K^ 

4 X^ cos('^~ + n9)K2 + 2X cos(“^ + nG) 

, 4 n cos(^ 4 n©)}K 2 ] 


(3.3c) 
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Equations (3.3) are rewritten in the following form 


by collecting the coefficients of U , V and W : 

n ’ n n 

00 


2 

n=1 

^n^12 ''^n®13 " ^”^12 



+ sin(^ 4- n9) = 0 

(3.4a) 

00 

2 

n=1 

^UnSl ^n"22 + ^^23 " ^/''n'^21 ^n'^22 




(3.4b) 

and, 



00 

2 

n=1 

^^n®31 ■'■ ^n®32 ■*■ ’'^'n®33 " ^ ^^n ^31 "'■ ^n^32 



+ W b 1 ] cosi-— + n9) = 0 

n oo 3 

(3.4c) 

where. 


®11 

= (A^, + B,\|) ^\6.* ®65 ' °66>"^ 

(3.5a) 

®12 

= + 2B^g + + (A^2 + ^66 ^ 2 


■f 

(3.5b) 

^13 

= (B,, + D,^)x2 + (3B,^ + D,^)a3i +[(B,2+2B,^ 

- 

DO^ 


+ A,2] X + (Bjg - + (Ajj - 

( 3.5c ) 

®21 

“ ^^16 ■'■ 2B, 5 + g) X^ + {A^2 + ^66 + B-i 2 + ^66) 

n 


'^26^ 

( 3.5d ) 
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®22 

^^66 ^^26 ^®26 026^^“^ 



+ (1 + ^22 )n^ 

( 3-oe) 

®23 




■ +■ + 2D2^)n2 + + B ] X + B22n^ +-n 

(3.5f} 

"31 




[(Bi 2 + 2Bg^ - + A^ 2^ ^ + (^2. - D2.)n^ 

^ ^^26 “ ®26 ^26 

(3.5g) 

"32 

= (Bi6 20^6^^'' + (B^2 + 2B^^ + D^2 -*- 



+ [(382^^ + 2D2^)n2 + A25 + 82^] ^ + B22n^ + n 

(3.5h) 

00 

00 

= 4\^n + [{D^2 + 2n^j3)n2 + b^^] 2^^ 



+ (2n^D2^ + 282^ - D2^)2nX + 022^“"^ “ ^ 

+ 822^ 2n2 _ 1 ) + 1 

(3.5i) 

and 

b^ ^ = a2k 2 + 2nX Kg + n^K^ 

( 3.6a ) 


b|2 - G 

(3.6b) 


b|3 = XK^ 

(3.6c) 


"^21"^ ^:12 ; ; ■ 

( 3. 6d ) 



( 3.6e) 
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23 

— 2 XK^ + 

nK.j 

( 3. of ) 

31 

"" ^13 


(3.6g) 

32 

^ ^23 


( 3 . 6 h ) 

33 

= ^2 + 

2nX Kg + n^K^ 

( 3. 6i J 


From equations (3.4), (3.5) and (3.6) it is evident 
that there are three equations ' corresponding to every integer n 
in the form: 


1 

®12 

^13 


u 

n 



*^12 

3 



®21 

^22 

^23 


V 

n 

- q 

“^21 

b 22 

^23 


V 

n 

^31 

®32 

®33 


W 

n 


^31 

*^32 

^33 


W ■ 

n 

— 


— 


— — 1 


— 


- ■ - 


; ^ 


(3.7) 

For numerical efficiency, and are eliminated 
from the above set of three equations, to get equation in 
Rewriting, the first two equations as; 


^12^n + ®224 + ®23'''n ‘ 



= 0 (3.8a) 
= 0 (3.3b) 


From equations (3.8), one has 

'/• - / . . - J— -j; {( 3^2823 “ ®22^13^n ’’^n 

^^1^22 - 

q [(322*^11 " ® 12 '^ 12 ^n ^'^n ^®22^12 " ^12^22^n \ 

+ ( 922^13 - ^ 12 ^ 23^0 %^ ^: ^^ 


(3.9a) 
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and 


V 


n 


(a a “■ ^23^11 ’'^n 

1^22 ~ ^I2^n 


q ~ ®12^11^n ^ ^^11^22 ®12^12^n '^n 


■'■ ^^11^23 " ®12’^13^n 


(3.9b) 


Equations (3.9) can be used for an iteration which aims 
at obtaining and in terms of W^. Iteration can be started 
by dropping the q terms, giving. 


U = ^^12^23 ^13^22^n 

/ ' ^^11^22 -^I2>n 


(3.10a) 


and, 


.. . ^®12 13 " ®11 23^n u, 

2~; n 

^®11®22 " ^12''n 


(3.10b) 


Substituting equations (3.10) in (3.9) we obtain and in 


terms of as: 


U 


1 


n 


^^11^22 " ®12^n 


^^^12^23 " ®22®13^n \ 
u \ ^ ^12^23 ~ ^13^22^n ^ 

q^ o 'v-i — — — — — ^ 


q [(322^11 - ai 2 ^ 1 2 ^n / -a^J 

^11 22 12‘^n 

/ K K 1 (a^2^l3 “ ^11^23.^ 

+ '®22 12 - h2 22>n ^^ 2 ^ 


n 


(a22'^13 - ®l2^23^n ’“^^n ^ ^ 


(3.11a) 


and 
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With these definitions, equations (3.11 ) are rewritten 


as : 


U 


n 


[ 


^1n r ^3n^1n ^4n^2n ^5n ^1 ,.t 


(3.12a) 


n 


n 


n 


and 


^n = 


r^2n r ^6n^1n ^n ^2n %n i •,■,» 

[— + + "2 + — ^ 


n 


R" 


n 


n 


(3.12b) 


n n 

Substituting equations (3.12) in the last equation of (3.7) and 
neglecting the higher order terms of q, the following form is 
obtained: 


R 


r 1 3 . f 4 in . 4n 2n , on ,-1 »i 

a. + q I — o — + o — + D — U 


^3^1 n , ^4n^2n 


bn 


‘13n 


n 


R 


R" 


R n 

n 


R 


2 n 


+ ^23n U, ^ 


n . ' n 

^6n^1n . ^n-^2n , ^n 


n 




n 


r 2 

n 


+ ^33n^n - ^ ‘’32n "“n + '^33n^’'’'n 


= 0 (3.13) 


This can be rewritten as 


[ ^3^„ , ^ ,,. i23 !2 , 3^^] . p[(b3 ,R, > b32R2 - >^33 


33"^ n 
n 


R^ - Rq)/R - (R3R^ + R 4R2 + ^6^ 


0 


(3.14) 


Equation (3.14) represents a set of n linear homogeneous 
equations in with q as the eigenvalue. These can be written 


as : 


[A ] IW! - q [b ] CWI . = 0 


(3.45) 
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= 


(a a - )~ “ ®23®11^n 'V 

^^11^22 ®12-’n 


( 3 ^ ^3 f 


a r(a ^ - h ^ '“12°23 " ^13®22^n ,.. 

q LU^^b ^2 - ® 12 '"l 1 '>n , 

^^1^22-^2^n 


n 


-^u\ i^io II 210 11 ur 

+ U^-,b22 “ ®12°12-'“ ■ 


^^12^13 ~ ^11^23^n 
b I ^ ^ ^2 ^ n 

^1 1 22 " 12^n 


+ ~ ^1 2®1 3^n.'‘'^n^ ^ 


(3.11b) 


Since further iteration of U and V^ yields an additional 

n n 

2 

expression in terms of q , which is negligible as treated by 
Ugural and Cheng [lb], equations (3.11 ) are taken as the final 
form. 

For computational convenience, the following quantities 
are defined as: 


and 


hn 


^^12^23 “ 

®22^13^n 

^2n 

=: 

^^12^13 ” 

®11®23^n 

^3n 

= 

^®22^11 " 

^12*^12^0 

^4n 

= 

^®22^12 “ 

^12^22^n 

%n 

=. 

^®22^13 “ 

^12^23^n 

%n 

'' = 

^®11^12 " 

®1 2^1 1 ^n 

^n 

■ =' 

^®1 1*^22 " 

®12’^12^n 

%n 


^®1 1*^23 " 

®12^13\ 



^®1 1^22 ” 

®1 2^1 2^n 


Subscript n in the above expressions implies the evaluation for 
specific value of n. 



where, [A] and [B] are diagonal matrices with elements and 
. given as : 


A = [■ 

y-. L 


n 


^13 ^1 ^23^2 , ^ -I 

/rsO -1 


R 


+ a 


33 


(3.16a) 


n 


B 


n 


r ^31^1 ^32^2 ^33 ^ ~ ^ 

" R 

r^sRi + ^4^2 ^^2 ] 


(3.1 6b j 


R- 


n 


and {WTf represents the column vector with elements .*• 


3.4 Optimization 

The objective of the present optimization study is to 
maximize the buckling strength of the laminated composite shell 
of specified geometry and weight. Such an optimization signi- 
fies the efficient use of composite laminates with different ply- 
angle orientations and ply thickness. The design variables here 
are fibre orientation and thickness of individual laminae. The 
number of layers of the shell is treated as a parameter, since, 
as a design variable it represents a mixed integer programming 
problem. Investigations are carried out for different loading 
combinations for the shell. 

3.5 Problem Formulation 

The optimization problem can be formulated as 

Maximize (f(o^^» 

) • ;: r- 

Minimize (-f(‘^j_» i — .1 , N 


where. 
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f(a,, h, ) : 


Eiih 


"til 

a. = orientation of i^ lamina 

*thi 

= thickness of i^ lamina 
N = No. of laminae 



Subject to the following constraints: 

(i) lower bound and upper bound on the fibre orientation 
angle 

0 ^ < 1 900 


(ii ) lower bound and upper bound on the lamina thickness h 

h. < h. < h. 

1 . - 1 “ 
min max 

Here q is the lowest eigenvalue obtained from the solution of 
C -P 

equation (3.15). 

3.6 Method of Optimization 

Since the objective function for laminated composite 
cylindrical shells is quite a complicated function of design 
variables a and h, constrained optimization method without using 
derivatives is suitable for solving the problem. The problem is 
to maximize the buckling load and find the fibre direction and 
lamina thickness at which buckling load is maximum. 

For the constrained optimization solution without using 
derivatives, sequential augmented La gran gi an multiplier method 
is used. This method finds the minimum of a multivariable non- 
linear function subject to non-linear bounds, equality and 
inequality constraints. The formal statement of the method 

can be given as [23] 
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Minimize f (X^ , 

X2, X3, . . x^) 


subject to 



\ < Gk ^^1 » ^2’ 

**•» ^n^ < ''^u 

k = 1 , M 

Gj^ (X^, X2, X3, 

• • • > ) > 0 • 0 

k = M+1 , MrfN 

Gj^ (X^ , X^, X3, 

0 

• 

0 

11 

c 

X 

• 

• 

• 

k = l'i+N+1 , A'k-N+P 


Optimization analysis was carried out using the programme 
E04UAF, supplied by Numerical Algorithm Group (NAG), available 
in the system library of DEC-1090 computer at IIT Kanpur. 

The bounds on the constraints adopted in our case are 
as follows: 

1. The angle of each ply is constrained between 0© and 90° 

2. The thickness of each laminae is constrained between 
0.01 in. to 0.1 in., which being the practical range of 
commonly available laminae thicknesses. 

3.7 Programming Procedure : 

The main steps of programming procedure to determine 
the eigenvalue are as follows: 

* For specified composite material properties and shell 

geometry the coefficients of stiffness matrices A^j, 
and are computed. 

* The coefficients aj^j and b^j given by equations (3.5) ano 
{3,6) are computed for each value of n ranging from 1 to 20. 

* Parameters a, defined by equations (3.12) are 

computed and the diagonal matrices [a 1 and [Bj given by 
equations (3. 16) are obtained^ 
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> Using NAG subroutine F02BJF available with DEC-1090 system 

the lowest eigenvalue solution of equation (3.15) is obtained. 
^ The procedure is repeated for values of m varying from 2 to 
20, and the lowest value of buckling load is determined. 

S'T'a Steps of Optimization Procedure 

The initial design values of and (i = 1 , N) are chosen 
within specified bounds and. the value of using OBJECT 
subroutine is calculated. 

^ For obtaining the maximum value of optimization is 

carried out using NAG subroutine E04UAF. 

The subroutine AMONIT is supplied for monitoring the optimi- 
zation procedure, and is called at the end of every cycle 
of E04UAF. 

^ Values of a^'s and h^^ ’ s at which the maximum value of q is 
obtained are given as the output along with the maximum 
value of q. 



CHAPTER 4 


M SULTS AND PISnn.R.’^TnM ' 

In this chapter, the results on the buckling and 
optimization of shells are presented and discussed. Eigen- 
value solution Of equation (3.15) is obtained to get the 
buckling load for a shell of specified geometry subjected to, 
all possible combinations of loadings. From this, a conclu- 
sion regarding the values of m and n, at which lowest eigen- 
value occurs, is made. Next the effect of laminate configu- 
rations on buckling loads is studied. For shells of specified 
geometry six different configurations are studied and the 
results compared with those by Khot [ll] who has used Donnell' 
•theory. The effect of number of shell laminae on buckling 
loads is studied for a shell of specified geometry subjected 
to axial and pressure loadings and compared with Jones and 
Morgan [14] results who used Donnell's theory. Finally the 
effect of ply-angle orientations on buckling loads is done 
for shells of specified geometry. Results are compared with 
the ones obtained by Whitney and Sun [16], who also used 
Cheng-Ho theory with a different solution technique. 

The computations have been carried out on DEC-1090 
computer system at I.I.T. Kanpur. The mate’rial properties 
of composite materials which are used for numerical computa- 
tions are given as follows: 
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Composite material 

X* 

; 

1 psi 
-1 

“T — 

: ^T 

1 

1 psi 

1 

i “^LT 

\ psi 

1 

^ 

1 I 

! LT ! 

} 1 

f 1 

! t 

TL 

Boron Epoxy 

30.0x10° 

3.0x10^ 

1 . 0 x 10 ^ 

CO 

• 

O 

0.03 

Graphite Epoxy 

30.0x10° 

0.75x10^ 

3.75x10° 

0.25 

0.00625 

Carbon Epoxy 

20 . 0 x 10 ° 

1 . 0 x 10 ^ 

0 . 6 x 1 0 ° 

0.25 

0.0125 


The definition of parameters and used in the 


tables is reproduced for convenience; 


= qK 


= 


and 


12 ~ '-i*'2 q^ ~ '^*'3 


= qK. 


All the buckling loads are obtained in terms of dimen- 


sionless 

parameter (qL /Ej^jt ) 

X 10' 

■, where 

q , 

= N , 


0 

and 

s = s = 

q 

= P , 

Ki 

0 

and 

II 

CO 

I! 

<N 

q 

= T , 

K 3 ^ 

0 

and 

li 

<N 

II 


4 . 1 Critical Buckling Load and Buckling Pattern ; 

Buckling loads are obtained for six layered Boron-Epoxy 
composite laminated shells under all possible loading combina- 
tions. The geometric properties used for this study are. 

shell ^ ^ ^ ^ 

Radius of shell * 

Thickness of each lamina ; 0.12 in 

Table 4.1 gives the value of non-dimensional buckling 
load for the shell for specified values of m and n under 
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pressure loading. For each value of m and n was varied from 
2 to 20 and m was varied from 1 to 20. The lowest buckling 
load of 100.84 x 10 is obtained for m = 18 and n = 12. In 
Figure 4.1 the variation of buckling load obtained as eigen- 
value of equation (3.15) for specified values of m and n, 
with changing values of m and n is -shown. Evidently, the range 

of variation of m and n chosen is sufficient to determine the 
buckling load. 

Tables 4.2 and 4.3 give similar results to determine 

the critical values of the non-dimensional buckling load for 

the shell under axial load and torsion respectively. For the 

shell under axial compression the value of critical buckling 

2 

load is 18.95 x 10 , occurring at m = 18 and n = 12. In the 

case of torsional loading, the value of critical buckling 

2 

load is 9,38 x 10 , occurring at m = 18 and n = 12. In Figures 
4.1 and 4.2 the buckling loads are plotted for different value 
of m and n. 

Table 4.4 gives critical values of non-dimensional 

buckling load for shell subjected to combined loading (K^ = 1, 

^2 = 1 » = 1 ). The non-dimensional critical buckling load 

2 

is found to be 2.91 x 10 at m = 18 and n = 10. In Figure 4.4 

the buckling loads are plotted for each value of ra and n. ^ ^ ^ 

computations are carried out for all possible 
cases of loadings and these result are concisely given in 
Table 4.5, From these results it is seen that for case of 
single loadings the buckling pattern is m = 18 and n = 1 2, 
and for case of combined loadings the buckling pattern is 
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m = 18 ano n = 10. This observation is very significant. It 
was later used in optimization problem which requires repeated 
determination of the eigenvalue. Using the values of m and n 
obtained here, direct comparisons are made with the results 
obtained by Khot [ll], discussed later. 

From Table 4.5 it is observed that torsional load 
affects the buckling of the shell most strongly. The effects 
of axial load and pressure load are in descending order. 

4 . 2 Effect of Laminate Configurations on Buckling Load : 


The shell considered for this' study has the following 
properties: 


Composite material 
Number of layers 
Length of shell 
Radius of shell 
Thickness of laminae 


Boron Epoxy 
3 

72.0 in 

6.0 in 
0.012 in. 


Types of configurations studied are: 
a) - a, 0, a b ) - a , Tt/ 2 , a c) 0, a, - a 

d) Tt/2, a, - a e) a , - a, 0 f) a, - a, Tt/2. 

Type of loading considered in all the cases is combined axial 

load, pressure and torsion. 

tables 4.6, 4.7 and 4.8 show the effect of the various 
configurations on the values of buckling loads. The results 
were obtained using m = 18 and n = 10. These results are 
compared with the results by Khot [1l] who used Donnell's 
theory of shells. 
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The results by Khot are 20% higher than the experi- 
mental ones. The fact that present results are more close to 
experimental results by about 1 0-1 2% implies that Cheng-Ho 
theory gives result which are significantly better than those 
obtained by Donnell's theory. Figures 4.5 and 4.6 show the 
variation of buckling loads with ply orientations for differ- 
ent configurations. The configurations - a, 0, a and - a, %/2,ci 
give maximum buckling load as. obtained by Khot at a = 60° and 
a = 300 respectively. This gives us an added confidence in 
the validity of Cheng-Ho theory and the inprovements it shows 
over Donnell' s theory. 


4.3 Effect of Number of Shell Laminae on Buckling Loads ; 

The shell considered for this study has the following 
material and geometric properties: 


Composite material : 
Length of shell : 
Radius of shell J 
Thickness of laminae : 
Orientation of laminae t 


Boron Epoxy 
34.64 in 
10.0 in 
0.12 in 
[0/90] g 


Table 4.9 gives the variation of buckling loads with 
the number of shell laminae. Two cases of single loadings - 
axial compression and pressure - are considered, i-igure 4.7 
shows the variation of non-dimensional buckling load with the 
number of laminae. It is observed that for pressure loading, 
the buckling load constantly increases with increase in number 
of layers; while for axial loading at N = 10, the buckling load 
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is msKimum sfisr wnich ii decrsssss. For diroctly coropsx'ing 
the present results with the ones obtained by Jones and J.'organ 
[ 14 ] the non-dimensional parameter for axial load is reduced to 
the form used by them. It is seen that the results by the 
present theory are less by about 15?^ as compared to the ones 
by Jones and Morgan who used Donnell's theory. 

4 • 4 Effect of Ply-angle Orientations on Buckling Loads : 


The shell considered for this study has the following 
material and geometric properties: 


Composite material 
Number of laminae 
Length of shell 
Radius of shell 
Thickness of laminae 


Carbon Epoxy 
4 

60.0 in 

6.0 in 
0, 1 2 in 


Table 4.10 gives the values of non-dimensional buckling 
loads for shells subjected to axial and torsional load singly. 
For the various configurations of ply-angle orientations, a 
reasonable agreement is found with the values of buckling loads 
obtained by Whitney and Sun [ 16 ] who used Cheng and Ho equations 
and obtained the buckling load by employing double Fourier 
series expansion for the displacement functions. The material 
and geometric properties considered here are identical xo the 
ones used by Whitney and Sun. The difference in the results by 
Whitney and Sun- and present theory show that the deviaxion or 
one resu.lts is significant on the lower side thus approaching 
the experimental conditions. Also the solution technique used 
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in the present work is a , single series expansion in which the 
size of matrices for computation -of eigenvalues is considerably 
reduced. This reduces computation time considerably. 

From the results of Sections 4.1, 4.2, 4.3 and 4.4 
the improvement of Cheng-Ho theory over Donnell’s theory is 
established for the ease of long shells. Comparison with the 
results obtained by Whitney and Sun validates the method used 
to obtain the solution using Gheng-Ho theory. Reasonable agree- 
ment with the available results [l6] justifies the use of m = 18 
and n = 12 for single loadings and m = 18, n = 10 for combined 
loadings.- 

Optimizatio n Re s ults ; 

Results are compared with the one' s obtained by Hirano 
[19] who has used Donnell's theory. Study is done for all 
cases of combined loadings, and for varying length- to- radius 
ratios of the shell. The effects of fixed lamina thickness and 
variable lamina thickness on maximum buckling load are also 
studied. 

Optimization study is done to maximize the buckling 
strength of laminated composite shells of speciiied geometry 
and weight. The design variables there are fibre orientations 
and ply thickness of individual laminae.. Number of layers of 
the shell is treated as a parameter. 

4.5 Optimum Buckling Loads for F ixed and Variable La . na jia 
C Thickness t 

The shell considered for this stody has the following 
material and geometric properties: 
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L - 60.0 in, L/a = 10.0, t = 0.8 in, N = 6. 

Table 4.11 gives the values of optimum buckling load 
for different length-to-radius ratios of the shell subjected to 
all cases of combined loadings. Each lamina is considered to 
be of equal thickness and the values of ply— orientation angles 
are obtained at the maximum buckling load. By varying length- 
to-raoius ratios of the shell it is observed that as this ratio 
increases the maximum buckling load decreases. This is due to 
increase in effective bending stiffness of shell. The values 
of buckling loads for the combined loading conditions yields 
the values of a which are less than those for loads acting 
singly. The general trend is similar to that observed in 
Section 4.1. 

Table 4.12 gives the values of optimum buckling loads 
for the shells with thickness of each lamina variable. At the 
optimal point the laminate configuration turns out to be very 

close to antisymmetric laminate. 

It is observed that the ratio t^/t^ ~ 3.05 and 
t./t^ 1.49 for each case of loading is independent of the 

length-to-radius ratio of the shell. For six layered laminates 
the optimal thicknesses ratios, satisfy the relation 

, Q_ A 
~ (R^ - 1 ) + f 1 ) 

where R 2 

for antisymmetric angle ply laminates. This was obtained by 
Shaima [23j to. predict decoupling in antisymmetric laminares. 
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4.6 Efle c t. °f.^tarting,_ Point on Optimal Ruckling Lnari.c; : 

The shell considered for this study has the fdllomng 
material and geometric properties: 

i 

Composite material : Boron-Epoxy 

N = 6, L = 60.0 in, L/a = 10.0, t = 0.8 in 

The effect of starting point on the optimum solution of 
a six— layered laminated shell for L/a ratio of 10.0 is shown in 
Table 4.13. The case of axial loading is considered. It is 
seen that with different starting points the fibre orientations 
of plies come out different, but the change in non-dimensional 
buckling load is very small. In the case of starting point 
s = 450 , the value of the load differs significantly. For 
comparing the results obtained by Hirano [19] who has used 
Donnell's theory, the non-dimensional buckling load parameter 
is reduced to the form used by Hirano [l9]. The results obtained 
by Hirano are 60-65% more than the experimental results. The 
present results show a significant improvement on Hirano' s 
results i.e. results are closer to experimental result by 15.7;a. 
Moreover Hirano has assumed the thickness of plies to be pre- 
assigned and equal, and found maximum buckling load, treating 
fibre orientation as the only variable. In the present worx, 
as given in Table 4. 14 thicknesses are also treated as design 
variables to search for an improved design. It is seen that at 
the optimal solution, laminate tends to become antisymmetric 

with thickness ratios to decouple the laminate. Hardly any 
effect of starting points on the optimal load is observed. Thick 
ness variation increases the load , carrying capacity by 10 127 q. 



Table 4.1 


98878 ’ 

Values of. critical buckling loads and buckling patterns for 
six-layered Boron-Epoxy composite laminated shell under 
pressure loading 


t 1 

t r 

m J n { 

1 t 

f f 

t t 

9 3 9 

qL VE^jt X 1 0 

p- 

t 

1 

f 

! 

f 

I 

1 

m 1 

t 

t 

j- 

1 

( 

n j 

f 

t 

X 

9 9 9 

qL /E^^t'^ X 10 

10 6 

111.83 


16 

6 

107.00 

8 

110.96 



8 

105.38 

10 

109.22 



10 

104.71 

12 

4 07.68 



12 

103.42 

14 

106.75 



14 

102.52 

16 

105.40 



16 

101.41 

18 

104.86 



18 

101 .21 

20 

104.42 



20 

101 .08 



— 

— - 





6 

110.43 

8 

108.01 

10 

107.14 

12 

105.53 

14 

104.82 

16 

1 04. 21 

18 

103.99 

20 

1 03.98 


6 

106.31 

8 

104.60 

10 

102.79 

12 

100.84 

14 

101 .47 

16 

101.41 

18 

102.87 

20 

104.82 


6 

1 08.57 

8 

1 06.54 

10 

1 06 . 1 3 

12 

104.97 

14 

103.93 

16 

102.51 

18 

102.12 

20 

1 01.98 


20 6 
8 

10 

12 

14 

16 

18 

20 


105.81 
103.27 
1 02.09 
101.86 
101 .01 
102.38 
103.71 
1 03.66 


Table 4.2 


Values of critical buckling loads and buckling pattern for 
six-layered Boron-Epoxy composite laminated shells under 
axial compression 



Table 4.3 


Values of critical buckling loads and buckling patterns for 
six-layered Boron-Epoxy composite laminated shells under 
torsion 



Table 4.4 


D /1 

**T 


Values of critical buckling loads and buckling pattern for 
six^layered^ Boron-.Epoxy composite laminated shells under 
combined axial, pressure and torsional loadings 
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Table 4.5 - 


-Critical buckling loads and buckling 
Boron-Epoxy composite laminated shell 
type of loadings 


pa'tterns for a six-layered 
subjected to different 





m 

J 

r 

n ! 

■ . t 

j 2 , ,3 2 

qL /Ej^-pt X 10 

1 

0 

0 

18 

12 

100.844 

0 

1 

0 

18 

12 

18.955 

0 

0 

1 

18 

12 

9.382 

1 

■ 1 

0 

18 

10 

18.767 

1 

0 

1 

18 

10 

19.103 

0 

1 

1 

18 

10 

6.778 

1 

1 

1 

18 

10 

2.914 
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Table 4.6 


Values of critical buckling load for three-layered Boron-Eooxy 
shell (K^ = 1, = 1, = 1 ) 


Configura tion 

r 

a 

J. 

Critical buckling load J 

(Ocr Ib/iri) I 

Percentage 

difference 

_ 

. . 


Results by J 

Khot [ll] 1 
! 

Results by | 

present study 1 

- _ . 1 


-a, 0, a 

10 

937 ■ 

• 

859.97 

8.27 

20 

1 062 

954.55 

10.16 

30 

1 000 

905.70 

9.43 

40 

1 000 

905,30 

9.47 

50 

1100 

978.78 

11.02 

60 

1 375 

1214.67 

1 1 .66 

70 

1250 

1 095 . 25 

12.38 

80 

1000 

896.73- 

10.33 

90 

900 

727.64 

0.04 

- a, 'n:/2, « 

0 

850 

747.06 

12.11 

10 

1000 

876.41 

1 2. 36 

20 

1200 

1057.08 

11.91 

30 

1300 

1138.02 

. 12.46 

40 

1150 

994.29 

1 3.54 

50 

1 050 

913.39 

13.01 

60 

1050 

916.33 

12.73 

70 

1150 

1033.16 

10.16 

80 

1050 

947.62 

9.75 

90 

lOOO 

928.80 

7.12 
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Table 4.7 

Values of critical buckling load for three-layered Boron-Eooxy 
shell = 1 , = 1, K3 = 1 ) ‘ 


Configuration 

Hf 

f 

t a 

f 

I 

t 

1 

! 

1 

-1 - 

- — ^ 

Critical buckling load ! 

Ib/in) | 

' T "" ■' ^ 

Results by| Results by J 

Khot [ 11 ] 1 present study; 

f I 

Percentage 

difference 

I 

•s 

0 

10 

512 ■ 

464.08 

9.36 


20 

567 

498.99 

11 .21 


30 

512 

463.72 

9.43 


40 

487 

448.72 

7.86 


50 

487 

450.62 

7.47 


60 

512 

461 .67 

9.83 


70 

567 

503.66 

11.17 


80 

512 

463.51 

9.47 


90 

■ 500 

470.25 

5.95 

7c/2, a, - a 

0 

500 

439.21 

12.16 


10 

515 

442.85 

14.01 


20 

575 

507.20 

11.79 


30 

515 

462.93 

10.11 


40 

475 

430.92 

9.28 


50 

475 

430.68 

9.33 


60 

515 

470.65 

8.61 


70 

• 575 

514.79 

10.47 


80 

515 

469.83 

8.77 


90 


454.72 

9.06 


% 
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Table 4.8 

-Values of critical budding load for three-layered Boron-Epoxy 
shell = "1 » ^2 ~ » ^3 “ ^ ^ 


Configuration 


Critical buckling load 

(Qcr Ib/in) 


Results by 
Khot [11J 


Results by 
present study 


Percentage 

difference 


a, - cc , 0 


426.74 

457.93 

451 .29 

463.78 

474.83 

460.43 

460.40 

455.15 

427.16 


10,16 
8.42 
11.51 
1 1 .86 
11.66 
9.72 
7.92 
8.05 
10.07 


- a, 11/2 


445.10 

437.85 

442.95 

457.94 
469.62 

473.96 
478 . 39 
459.79 
460.76 

453.95 


7.28 
9.73 
11.41 
11.08 
12.22 
1 2.23 
10.58 
10.72 
8.76 
9.21 
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Table 4.9 


Variation of critical buckling load with number of shell 
laminae for cross-ply laminated Boron-Epoxy shells 


Number of 
layers 

- 1 


f 

9 ’ 

X 10 ! 

f 

- 1 

Percentage 

difference 

2 ' 3 

pL X 10“^ 

- 

t 

Results by { 
Jones and ! 
fbrgan ! 

[14] I 

t 

Present I 
results . j 

I 

t 

f 

I 

t 

3 

17.33 

14.73 

14.98 

94.16 

5 

18.65 

15.85 

15.1 

98.73 

7 

20.15 

17.25 

14.37 

1 07 . 22 

9 

20.38 

17.87 

12.29 

116.70 

10 

20.83 

17.85 

14.31 

118.39 

20 

20.18 

17.82 

11.73 

131 .44 

25 

20.08 

17.66, 

12.05 

■ 1 38.21 

30 

19.51 

^ t6.56 

15.11 

146.08 

50 

i , 

18.72 15.87 

'i;,. 

15.23 

172.86 


so 


Table 4*10 


Variation of buckling load with ply— angle orientations for four- 
layered carbon epoxy laminated shells 


Laminate 
confi gu- 
ration 


Mon-dimensional buckling load (qL^/Ej^^t^) x 10 ^ 

- ___ 

Axial load J Torsional load 


t ^ " '' ■' f ' ^ ■ ■■ - - r , - 

{Results {Present {Percentage {Results {Present {Percentag 

{by {results {difference {by {results {differenc 

r^Vhitney 8 . { { ■ {Whitney & { ' 

ISun [16] { { {Sun [16] { { 

-J 1 1 t 1 


[ 0 / 0 ] 3 

11.6 

12.0 

+3.5 

1 .8 

1 .6 

-1 2.5 

[ 1 V 0 ] 3 

15.0 

15.2 

* +1 .3 • 

1 .9 

2.0 

+ D . 0 

[30/0] g 

18.4 

17.4 

-5.75 

2.8 

3.2 

+ 1 2.9 

[ 45 / 0 ] 3 

22.2 

20.8 

-6.73 

6.0 

5.0 

- 20.0 

[60/0]g 

23.6 

21 .2 

-11.3 

7.1 

6.4 

+1 1 .3 

[ 75 / 0 ] 3 

21 .0 

20.5 

-2.4 

9.1 

8.0 

-13.75 

[90/0] 3 

18.0 

18.2 

+ 1.1 

8.4 

7.8 

-7.69 
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Critical buckling load (qL^/Elj l^)x 10 



Critical buckling load ( qL^/ELT*^^’* 





Critical buckling load ( qL /e^j t ) x 10 




Critical buckling load ( qL^/ELj t'’)x10 





Critical buckling load (qcr)lb/in 



FIG. 4 5 ^EFFECT OF VARIOUS LAMINATE CONFIGURATIONS ON CRITICAL 
BUCKUNG LOAD FOR 3- LAYERED BORON-EPOXY SHELL 



Critical buckling load (qcr^ Ib/in 


70 



FIG. 4 6 EFFECT OF VARIOUS LAMINATE CONFIGURATIONS ON CRITICAL 
BUCKLING LOAD FOR 3-LAYERED BORON-EPOXY SHELL 



FIG. 4.7 EFFECT OF NUMBER OF SHELL LAMINAE ON CRITICAL BUCKLING LOAD 
FOR CROSS- PLY LAMINATED BORON-EPOXY SHELL 






Non-dimensional buckling load t^)x 10 
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Angle of ply orientation ( «^ ) 

FIG. 4.8 EFFECT OF PLY- ANGLE ORIENTATIONS ON CRITICAL BUCKLING 
LOAD FOR 4 -LAYERED CARBON EPOXY LAMINATED SHELL 
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CHAPTER 5 
CONCLUSIONS 


On the basis of the results obtained on the present 

work, the following general conclusions can be drawn: 

* The buckling of six layered composite laminated shells 
subjected to single loadings is found to occur at values 
of m = 18 and n = 12. In case of combined loadings, 
buckling is found to occur at values of m = 18 and 

n = 10. The solution technique adopted does not allow 
the satisfaction of any boundary conditions. Thus 
Cheng and Ho theory with this form of solution tech- 
nique can be used for the analysis of laminated compo- 
site shells, with good reliability. 

* By studying the effect of various laminate configura- 
tions on buckling load and comparing the results by 
those obtained by Khot [ll] who has used Donnell's 
theory, it is concluded that the use of Cheng and tfo, 
theory gives results closer to experimental values by 
about 10-15^, as compared to the results obtained using 
Donnell's theory. 

* The study of effect of number of cross— ply shell laminae 
on buckling loads indicates that for the case of axial 
loading, with the increase of number of laminae upto 10, 
the buckling load increases, after which it falls. In case 
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of pressure loading the buckling load increases with 
increase in number of laminae. 

* The study of effect of ply-angle orientations on buckling 
loads indicate a reasonably good comparison of the present 
results with those by ’.Vhitney and Sun [l6] who have used 
Cheng-Ho theory with the solution in form of double Fourier 
series expansion. The major advantage of using the solu- 
tion technique adopted in the present work is that computer 
time for numerical calculations is greatly reduced, and 
the results obtained are closer to the experimental values. 
As length- to-radius ratio of the shell increases, the 
m*aximum buckling load of the antisymmetric six-layered 
shells decreases. 

* For six-layered antisymmetric angle-ply laminated shells 

subjected to all combinations of loadings, the thickness 
ratios « 3.05 and t^/t^ ~ 1.5 have constant values 

independent of length- to-radius ratio and loading condi- 
tions. kX the optimal point laminate is antisymmetric 
with thickness ratios producing a decoupled laminate. 

Thus an antisymmetric, uncoupled six-layered laminated 
shell has the maximum buckling load carrying capacity. 

* Effect of starting point on optimal solution indicates 

the existence of different sets of ply orientations leading 
to almost the same buckling load. 

* In case of laminated composite shells, best fibre direc- 
tions cannot be obtained, but it is seen that at the opti- 
mum point, the values of maximum buckling load are constant. 
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* Optimization of buckling load with ply-thickness and 

orientations as design parameters increases buckling load 
by about 12 % over the case when laminae thickness is 
uniform. 




76 


LIST OF REFERENCES 


1. Budiansky, B., ’’Notes on non-linear shell theory", Journal 
of Applied M,echanics; Vol. 35; 196S; pp. 393-402. 

2. Mah, G. B. , Almroth, B.O., Pittner, E.V., "Buckling of 
orthotropic cylinders”, AIM Journal; Vol. 6; 1968; on. 
598-602. 

3. Tennyson, R.C. , "Buckling of circular cylindrical shells 
in axial compression", AIAA Journal; Vol. 2 \ 1964; pp. 
1351-1353. 

4. Becker, H., Gerard, G., "Elastic stability of orthotropic 
shells", Journal of Aerospace Sciences; Vol. 29*, 1962; 
pp. 505-512. 

5. Dong, S.B, , Pister, K.S., Taylor, R, L., "On the theory of 
laminated anisotropic shells and plates". Journal of 
Aerospace Sciences; Vol. 29; 1962; pp. 969-975. 

6. Tsai, J., "Effect of heterogeneity on the stability of 
composite cylindrical shells under axial compression", 

AIAA Journal; Vol. 4; 1966; pp. 1058-1062. 

7. Serapico, J.C., "Elastic stability of orthotropic cylind- 
rical shells subjected to axisymmetric loading", AIAA 
Journal; Vol. 1 ; 1963; pp. 128-138. 

8. Lowe, G.G. , "Structural analysis of orthotropic shells", 
AIAA Journal; Vol. 4; 1966; pp. 1843-1849. 

9. Dow, N.F., Rosen, B.W. , "Structural efficiency of ortho- 
tropic cylindrical shells subjected to axial compression", 
AIAA Journal; Vol, 4; 1963; pp, 481-490, 

10. Weingarten, V.I., "Effects of internal pressure on the 
buckling of circular cylindrical shells under bending". 
Journal of Aerospace Sciences; Vol. 29; 1962; pp. 804-813. 

11. Khot, N.S. , "Buckling behavior of composite cylindrical 
shells under axial compression", AIAA Journal; Vol. 8; 
1970; pp. 229-239. 

12. Cheng, S., Ho, B.P.C., "Stability of heterogeneous aleo- 
tropic cylindrical shells under combined loading", AIAA 
Journal; Vol. 1; 1963; pp. 892-898. 

13. Cheng, S., Ho, B.P.C., "Some problems in stability of 
heterogeneous aleotropic cylindrical shells under combined 
loading", AIAA Journal; Vol. 1; 1963; pp. 1 603-1 60o. 



77 


14. Jones, R.M. , fvbrgan, H.S., "Buckling of cross-ply lamin- 
ated circular cylindrical shells", AIAA Journal; Vol. 13; 
1975; pp. 664-673. 

15. Ugural, A.C. , Cheng, S. , "Buckling of composite cylindrical 
shells under oure bending", AIAA Journal; Vol. 6; 1968; 

pp. 349-354. 

16. ’Vhitney, J.M. , Sun, C.T., "A refined theory for laminated 
anisotropic cylindrical shells", Journal of Composite 
Materials; Vol. 41 ; 1975; pp. 471-482. 

17. Nshanian, Y.S., Pappas, M., "Optimal laminated composite 
shells for buckling and vibrations", AIAA Journal; Vol. 21 ; 
1983; pp. 403-438. 

18. Hirano, Y., "Buckling of angle-ply laminated circular 
cylindrical shells", Journal of Applied Mechanics; Vol. 46; 
1980; pp. 233-234. 

19. Hirano, Y.. "Optimization of laminated composite plates 
and shells". Composite ?i(1aterials - Recent Advances; 

Pergamon Press, N.Y. * 

20- Fliigge, W. , "Stresses in shells", Springer-Verlag, Berlin; 
pp. 433-483. 

21. Calcote, L. R. , "The analysis of laminated composite struc- 
tures", Van Nostrand Reinhold Company, New York; pp. 143- 
151. 

22. Jones, R.M. , "Buckling of circular cylindrical shells with 
different modulii in tension and compression", AIA'^ Journal; 
Vol. 9; 1971; pp. 53-57. 

23. Sharma, S. et.al., "A study of coupling in lam.inated plates", 
Fibre Science & Technology, An International Journal; 

Vol. 18; No. 4; 1983. • 



A£ - 1 ‘j (3 ^ - fn~y^fY)- cpt 


■®*nS878 



